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Abstract 

Starting from correlation identities for the Blume-Capel spin 1 systems and using cor- 
relation inequalities, we obtain rigorous upper bounds for the critical temperature. The 
obtained results improve over effective field type results. 



1. Introduction 

Correlation inequalities combined with exact identities are useful in obtaining rigorous 
results in statistical mechanics. Among the various questions that are resolved by them 
one is the decay of the correlation functions. The decay of the correlation functions give 
information about the critical couplings of statistical mechanics models. In this work, the 
method will be applied to study systems described by the spin one Blume-Capel model 
0,0. Firstly, we present the derivation of an exact relation for the two spin correlation 
function, valid in any dimension, which is an extension of Callen's identity for spin 1/2 
Ising model Q. Starting from these identities we will then make use of the first and 
second Griffiths inequalities and Newman's inequalities to obtain the exponential decay 
of the two spin correlation function . The coupling constant which are the upper bounds 
for the critical temperature are obtained for d=2 and d=3 dimensions. In this study 
the coupling parameters obtained improve effective field results. Upper bounds for the 
critical temperature T c for Ising and multi-component spin systems have been obtained 
by showing (for T > T c ) the exponential decay of the two-point function 0, H 0. Spin 
correlation inequalities and their iteration are used by Brydges et al Q, Lieb i[7| and 
Simon The procedure to improve the bound for the critical temperature over the 
effective field result for the classical 5=1 model is as follows: starting from a two-point 
correlation function identity, a generalization of Callen's identity |3|| for this mod el [p 



and using Griffiths 1st and 2nd inequalities (Griffith I, II) (see [9J, [I0j.[llj.[12l|.|l3j.|l4j) 
and Newman's inequalities 1(J HI we establish the inequality for the two-point function, 
< SqSi >, as 

< S Si >< J2 a j < SjSi >, < a,j < 1 (1) 
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which when iterated (see [5j) implies exponential decay for T > T c . In section 2 we 
present the derivation of the correlation identities for the Blume-Capel model Q. In 
section 3, we apply these identities to the d = 2 and d = 3 lattices. Next, in section 4, we 
apply the correlation inequalities to obtain the upper bounds for T c . Numerical results 
can be found in section 5, and in section 6 we present our concluding remarks. 

We write the Hamiltonian for the classical spin one system, known as the Blume-Capel 
model, as 

H = -jJ2s i S j -DY / Sl (2) 

i,j i 

where J > 0, D is the single ion anisotropy and the first sum is over the nearest neighbours 
spins on the lattice. We define the thermal average < ... > by 

<... >=Z- 1 Y / (-)e- pH ,Z^Y, e ' PH ( 3 ) 

{St} {S z } 

where each Si is restricted by Si — — 1, 0, +1. 



2. Correlation identity for the spin one model 

We reproduce the generalization of Callen's identity for the spin 1 Blume-Capel model 
which has been obtained previously by Siqueira and Fittipaldi [8J. Let 

, prose ^_ Tr(F(S)S t e-^) 

where F(S) is any function of S different from We can write H — Hi + H', where 

ffi = -(EWi-^ (5) 

b'l=i 

is the Hamiltonian describing site i and its neighbours, and H' corresponds to the Hamil- 
tonian of the rest of the lattice. Consequently [Hi,H'] = 0. From Eq.@ and Eq.©, we 
get, 

, EVcw ^_ TrF(S)e-^ H -+ H '^S l _ Tr'TnF(S)e-^- S t e^ H ' 

Tre -p(H i+ H>) Tr>Tne-0 H 'e-0 H ' [) 

or 

< >= ^ Wr ^ (7) 

where Tr'Tri = Tr. Finally, we obtain, 

<F{S) Si >={F{S) lr * e _J; ). (8) 



Explicitly operating the trace Tr,-, we get 
<F{S)S l > = (F(S) 



2e^ D sinh{Y J] PJ l3 S j ) 



2eP D co S h(J2 i l3J lJ S J ) + l 



F(S) J] e^ s ^)f(x)\ x=0 , (9) 
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with V = -ff-, such that e f(x) — f(x + a), and 

,, . 2e /3D sinh(x) , , 

™ = 2e^cosh {x) + l ^ 

As Sf l = S] and Sf l+1 = Sj for n = 0, 1, 2, 3, we obtain, 

= S'|cosfr(A) + S>smfr(A) + l-S J 2 , (11) 

and, applying Eg. (ITU)) and Eq.JTJJ in Eq.©, we get 

< F(S)S l >= (F(S) J] {S 2 cosh{(iJ l3 V) + S ]S inh{pJ t] V) + l-S 2 ))f(x)\ x=Q (12) 

Similarly for the correlation function involving the square of the spin function Sf, we 
obtain, 

< G(S)Sf >= (G(S) J] e^ v )g(x)\ x=0 , (13) 

|j|=i 



with, 



resulting in, 



. , 2e fjD cosh(x) 



<G{S)S 2 >=(G{S) Yi {S 2 j cosh{pJ i:j V)+S j sinh{pJ ij V)+l-S'l)^g{x)\ x= Q. (15) 

37^,1.71=1 

The function G(S') is any function of S, except Sf . The equations (IT21 and (fT5|) are exact 
and generalize Callen's identity which was obtained for the S = 1/2 Ising model [3j. 

3. Exact correlation identities applied to the d=2 and d=3 lattices 

Let us apply the previous results for < F(S)Si > and < G(S)S 2 > given by equations 
(|12[) and (|15p for specific lattices in two- and three-dimensions. The two spins correlation 
functions, < SqSi >, are obtained from equations (1X2"]) and (ITS]) by defining i^(5) = S*;. 

For the d — 2 and z = 3, i/ie honeycomb lattice 
We obtain from Eq. (jT3 



+A 3 < S l S 3 S k S l > +A 4 < S l S 2 S 2 k S l >, (16) 

i<j<k i<j<k 
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where the A coefficients are given in appendix | Appendix A.l We also obtain, from 
Eq.CES, 

< S$S l > = B a + B 1 J2 < SfSi >+B 2 ^< S l S 3 S l > 

i i<j 

+B 3 J2 < S i S j S i > + B ± J2 < S i S i S kSi > 

i<j i<j<k 

+B 5 J2 < SfS^SlSi >, (17) 

i<j<k 

where the B coefficients are given in appendix | Appendix A.l| 

3.2. For the d — 2 and z — 4 7 the square lattice 

We obtain from Eq. (fT2")) for the two spin correlation functions < SqSi > the expres- 
sion, 

< SoSt > = Ai ^ < SiSi > +A 2 < SiSjSl >+ A 3 < S * S J S kSi > 

i i<j i<j<k 

+A4 < SiSjS^Si > +A§ < SiSjSkS^Si > 

i<j<k i<j<k<m 

+A 6 ^ < SiS?jSlSf n Si >, (18) 

i<j<k<m 

where the A coefficients are given in appendix |Appcndix A.2| We also obtain, for the 
function < SqSi >, 

< S%S t > = B + B 1 J2< SfSi >+B 2 ^< SiSjSi >+B 3 J2< s ? s ] s i > 

i i<j i<j 

+B4 y ' < SiSjSfrSi > +-B5 y ' < Sf Sj SfcSi > 

i<j<k i<j<k 

+B$ ^ < SiSiSkS m > +B 7 ^ < SiSjSlSf n > 

i<j<k<m i<j<k<m 

+B S ^2 < SfSjSlS^Si >, (19) 

i<j<k<m 

where the B coefficients are given in appendix |Appendix A.2 



3.3. For the d = 3 and z = 6 ; the cubic lattice 
We obtain from Ea. (fT2l) 



< SoSi > = At < StSi >+A 2 ^< SiS]Si >+A 3 < SiSjS k Si > 

i i<j i<j<k 

+A4 < SiSjSj:Si > +A$ < SiSjSkS^Si > 

i<j<k i<j<k<m 
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+Ag S ' < SiSj Sj c S m Si > -\-A-j S ' < SiSjSkS m S n 

i<j<fc<m i<Cj <k<m<n 

+Ag < SiSjSkS^S^Si > 

i<j</c<m<n 

> 

i<Cj </c<m<n<p 

H~^-io ^ ^ < SiSj Sf t S rn S n Si 

i<j <k<m<n 

+A n J2 <S l S :l S k SlS 2 n S^S l >, (20) 

i<J<fc<m<n<p 

where the A coefficients are given in appendix |Appendix A.3 We also obtain, for the 
function < S^S/ >, 



< Sfa > = B + B 1 < SfS t >+B 2 J2< S l S J S l >+B 3 J2< s f s ] s i > 

i i<j i<j 

+£?4 y < SiSjS^Si > +-B5 y < sf Sj SfcSi > 



i<j<k i<j<k 

+ i?6 < SiSiSkSm > +-S7 < SiSjS^Sm > 

i<j<fe<m i<j < fc<m 

i<j<k<m 

+ -~^9 < SiSjSkS m S n S p Si > 

i<j</c<m<n<p 

+-S10 ^ < SfSjSlS^ L S^Si > 

i<j <k<m<n 

~\~ B \ \ ^ ^ <c (S'j tS' ~~ S m S n Sp S[ ^> 

i<ij <fc<m<n<p 

+-B12 < SfSjSlS m S n Si > 

i<j <k<m<n 

+-^13 y ] < Sj S^S^SnSpSi > 

i<j<fc<m<n<p 

+ ^ 14 5^ < S?SjSkSmS n Sl > 

i<j <fc<m<7i 

~^^ 15 < $i Sj SkSmSnSpSl >, 

i<j<fc<m<n<p 



(21) 



where the B coefficients are given in appendix |Appcndix A.3 



The sums over i, j, k, m, n and p are over the nearest neighbors of to which we have 
given a numerical ordering. The proof of results (1161) for the case (|3.ip . the honeycomb 
lattice, is presented in the appendix |Appcndix B| as an example for the other cases. 



4. Application of the correlation inequalities 



In the following results we will made use of the following inequalities: < S A >> 
(Griffiths I), < S A S B > - < S A >< S B >> (Griffiths II) (see @, 0,0, [3), 



< SiF >< Y,j < S l S J >< dF/dSj > (Newman's) ([H, UdDand < SfS A ><< S A > 
([12|>[l3l), where < S A >= YliSi, < Sb >— Yli^i and F is a polynomial function of 
variables S. 

From the equations for the two spin correlation functions obtained in subsections 
(|3.ip , (|3.2p and (|3.3I) and applying the Griffith's and Newman's inequalities we obtain 
an inequality of the form 

< S S t >< < S ^ >' ( 22 ) 

1*1=1 

where cij is a sum of products of two-point functions. 

(a) Case d = 2, z = 3, honeycomb lattice. 
Using 

< S?SiSi ><< > (23) 
in equation (|19[) . in the A 2 term and Griffiths II , i.e., 

< SiSjSkSi >>< SiSj >< SkSi > (24) 

in the A3 term, and noticing that A 2 and A3 are negative, we get for d=2, z=3, 

< S Si >< {A x - I A 2 I - I A 3 |< >id +A 4 ) J2 < s i s i >, (25) 

|i=i| 

(b) Case d = 2, z = 4, square lattice. 

Using inequality (|2"3"]) in equation (|20p . in the A2 term (A2 < 0) term, Griffiths II in 
the A3 term (A3 < 0), the inequalities 

< SplSiS! ><< SiSi > (26) 

in the A4 term (A4 > 0) and 

< SplSiS^ ><< > (27) 
in the Aq term (A@ > 0) and in the A5 term using Griffiths II, we get for d=2, z=4, 

<S S l > < (A 1 -\A 2 \-<S 1 S 2>1D \A 3 \ 

+A 4 + < S 1 S 2 > 1D A 5 + A 6 ) ^ < S i S l >, (28) 

|<=i| 
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(c) Case d = 3, z = 6, cubic lattice. 



As before, we use in equation (|20p. inequality ((23"]) in the A2 term (A2 < 0) term, 
Griffiths II in the A3 term (A3 < 0), the inequalities ([2"Bf in the A 4 term (Aj > 0), 
inequality (24) in the A§ term (A > 0) , and Griffiths II in the A$ term. For the term 
A 7 (> 0) we use Newman's inequality and for the terms A$(> 0), Ag(> 0), Aiq(> 0) and 
An(> 0), we use inequality (22). Then, we get for d=3, z=6, 

<S S l > < (A 1 -\A 2 \-<S 1 S 2 > 1D \A 3 \+A 4 
+ < S1S2 >id A 5 

+A 6 + A 7 + A 8 + A 9 + A w + A n ) < SiSi > (29) 

|i=i| 

The two-spin correlation function < S1S2 >id is the one-dimension model two spin 
correlation separated by a distance of two lattice sites. By bounding the resulting two- 
point function occurring in the previous results from below with the two-point function 
of the one-dimensional infinite chain (see Appendix | Appendix B ) , we get: 



< S Si ><Y, a *< S * S i >. ( 3 °) 
|t=i| 

where, 

(a) For d = 2, z = 3, honeycomb lattice. 

aj =A 1 -\A 2 \-\A 3 \< SoSi > 1D +A 4 ; (31) 

(b) For d = 2, z = 4, square lattice. 

aj = Ai- \A 2 \-< S!S 2 >id\ A 3 I +A 4 

+ < SiS 2 >id A 5 + A 6 ; (32) 



(c)For d = 3, z = 6, cubic lattice. 

a = Ax- \A 2 \-< SiS 2 >id\ A 3 I +A 4 
+ < SxS 2 >xd A 5 + A 6 

+A 7 + A 8 + A 9 + A 10 + An • (33) 
The one-dimensional correlation function is given by (see Appendix |Appcndix C[ ) : 



l + y/(l- 2/(2/3 J) 

< SlS2 >1D - J(2pT) (34) 



and /(2/3J) is given by fll 
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Table 1: Estimatives for kTc/J for D = in previous and in the present work. 





d = 2,z = 3 


d = 2,z = 4 


d = 3, z = 6 


MFA 


2 


2.667 


4 


Siqueira 


1.518 


2.188 


3.516 


Yuksel 




1.964 




CVM 






2.886 


Series 




1.688 


3.192 


RG 




2.128 


3.474 


Monte Carlo [30] 




1.695 




Monte Carlo [29] 




1.681 




Wang Landau [28] 




1.714 




Present work 


1.591 


2.322 


3.678 



Table 2: Estimatives for kTc/J for D = oo in previous and in the present work. 





d=2, z = 3 


d = 2,z = 4 


d = 3, z = 6 


MFA 


3 


4 


6 


Siqueira 


2.103 


3.088 


5.076 


CVM 






3.876 


Series 






4.482 


RG 




2.884 


4.932 


Present work 


1.999 


3.070 


5.084 



5. Numerical results 

Evaluating numerically the value of T such that J2 a j — 1) a i > 0, we obtain, by 
sufficient condition, upper bounds for T c , which are shown in tables[T]and[2J m comparison 
with results obtained by other methods. 

For the evaluation of the self-correlation terms (< Sf >) that emerge from the applica- 
tion of the Griffith's and Newman's inequalities, we use, for the D = case, < Sf >< 2/3, 
correct for a spin 1 ferromagnetic system, and, for the D = oo case, < Sf >= 1, since in 
this limit the Si = spin value is suppressed. 

For the honeycomb lattice our result has to be compared with the mean field and 
the effective field calculations. Those results are not rigorous, as ours, and the numerical 
values we obtain improve those mean field type results and therefore represent the upper 
bounds, for the limits D = and D = oo. For the square and cubic lattices besides the 
mean field type results, for which the previous comments apply, there are other results, 
better than mean field type, obtained by series and renormalization group calculations, 
which can be used as a comparison. The importance of the present numerical results lies 
in the fact that they were obtained using an identity and rigorous inequalities for the 
two-spin correlation function. For this reason they represent rigorous upper bounds for 
the critical temperature. 
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6. Final Comments 



We have presented the derivation of correlation identities for the Blume-Capel spin-1 
model which are exact in all dimensions, and we have made use of correlation inequalities 
to obtain the upper bounds for the transition temperature. The coupling constants 
obtained for those bounds are calculated for d=2 (honeycomb and square lattices) and 
d=3 (cubic lattice). We obtain rigorous results that improve mean field type calculations. 

Appendix A. Coefficients of the Spin Correlation Identities for d=2, z—3 
and z—4. 

With k = j3J and f{x) given by relation (fTU|) . we have for 
Appendix A.l. d=2, z=3 

A 1 = 3f(k) > 0, (A.l) 
A 2 = {if {2k) - 6/(fc)) < 0, (A.2) 

A 3 = ~(f(3k)-3f(k)) <0 (A.3) 
M = \ (5/(fc) + f(3k) - if {2k)) > (A.4) 



and 



B =.9(0), (A.5) 

B 1 =3{g{k)-g{0)), (A.6) 

B a = |G7(2fc)-5(0)), (A.7) 

B 3 = ^g{2k) + -6g{k) + ^g{0), (A.8) 

B 4 = ^{g{3k) - g{k) - 2g{2k) + 2g{0)), (A.9) 

B 5 = \ 9 {ik) - *-g{2k) + ^g{k) - ~ 5 (0). (A.10) 



Appendix A.2. d=2, z=4 



A 1 =4f{k)>0, (A.ll) 

A 2 = 6/(2*;) - I2f{k) < 0, (A.12) 

^3 = /(3*0 - 3/(fc) < 0, (A.13) 

Ai = 15/(fc) - 12/(2*0 + 3/(3*)) > 0, (A.14) 

A 5 = if (4k) - f(3k) - }{2k) + 3f(k) > (A.15) 
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A 6 = i/(4A) - 3/(3*) + 7/(2*) - 7/(*) < (A.16) 

and 

B - 5 (0), (A.17) 

Bi=4( fl (fc) - 5 (0)), (A.18) 

B 2 = 3(5(2A:)- 5 (0)), (A.19) 

B 3 = 3( fl (2fc) - 4 5 (fc) + 3.g(0)), (A.20) 

B 4 = 3(5(3*) - 2 5 (2fc) - g(k) + 2.g(0)), (A.21) 

B 5 - g(3k) - 6g(2k) + 15g(k) - 10 5 (0), (A.22) 

B 6 = l -{g(Ak) - 4g(2fc) + 3.g(0)), (A.23) 

B 7 = |g(4*) - 3 fl (3fc) + 3.g(2fc) + 3g(k) - ™g(0), (A.24) 

17 35 

B 8 = -5(4fc) - ff (3fc) + -5(2*0 - 7g(k) + —g(0). (A.25) 

Appendix A. 3. d=3, z=6 

A\ = 6/(fc) > 0, (A.26) 

A 2 = -30/(*) + 15/(2*) < 0, (A.27) 

A 3 = 5/(3*) - 15/(fc) < 0, (A.28) 

At = 75/(fc) + 15/(3*) - 60/(2*) > 0, (A.29) 

A 5 = -15/(3fc) + 45/(fc) + y/(4fc) - 15/(2*) > 0, (A.30) 

A 6 = -45/(3*) - 105(/(*) - /(2*)) + y/(4fc) < 0, (A.31) 

3 15 15 

A 7 = g/(5fc) - -/(3fc) + T /(fc) > 0, (A.32) 

45 105 15 

A 8 = -f(3k) - — /(*) + -f(5k) - 15/(4*) + 30/(2*) < 0, (A.33) 

3 15 15 3 3 15 

A 9 = --/(5*) + -/(3*) - -f(k) + -/(6fc) - -/(4*) + -/(2fc) < 0, (A.34) 

405 315 15 

Ao - -jp/(3fc) + ^/(fc) + y/(5fc) - 15/(4*) - 90/(2*) > 0, (A.35) 

5 45 15 135 15 5 

A n = --/(3fc) + -f{k) + -f{Ak) —f(2k) -f(5k) + -/(6fc) > (A.36) 

and 

B = 5(0), (A.37) 

S 1 = 6(5(*)-5(0)), (A.38) 
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5 2 = _(.g(2fc)-.g(0)), (A.39) 

B 3 = H ( g (2k) - 4.g(fc) + 3.g(0)), (A.40) 

B 4 - 15(.g(3fc) - 2 5 (2fc) - g(fc) + 2g(0)), (A.41) 

5 5 = 5( S (3fc) - 6g(2fe) + 15.g(fc) - 10g(0)), (A.42) 

B 6 = H ( 9 ( 4 fc) - 4 ff (2fc) + 3.g(0)), (A.43) 

S 7 = 45(i 5 (4/c) - g(3k) + g(2k) + g(k) - ^(0)) , (A.44) 

1 7 35 

B 8 = 15(- 5 (4fc) - 3(3*0 + 2-9(2fc) - 7.g(fc) + y 5 (0)), (A.45) 

B 9 = ^(.9(6fc) - 6 5 (4fc) + 15g(2fc) - 10.g(0)), (A.46) 

B w = f (-1260(0) + 45g(3fc) + 210g(fe) - 120g(2fc) - 10.g(4fc) + g(5k)), (A.47) 
8 

3 1 fi^i 77 1 11 

= -(- Y g(3k)-66g(k) + —g(2k) + Y g(0) + -g(6k) + Y g(4k)-g(5k)), (A.48) 

#12 = ^(-8.g(2fc) + 140(0) + 0(5*;) - Ug(k) + 13g(3k) - 6g(4/c)), (A.49) 

Bi 3 = ^(-40.g(3fc) + 48 ff (fc) + 15.g(2fc) - 42 5 (0) + 26.g(4fc) + g(6k) - 8g(5k)), (A.50) 

B u = — (-20(4fc) + 80(2fe) - 60(0) + ,g(5fc) - 3.g(3fc) + 2.g(fc)), (A.51) 
8 

Bin = :|( 2 s( 4 *0 - 170(2fc) + 140(0) - 4g(5k) + I2g(3k) - 8g(k) + g(6k)). (A.52) 

Appendix B. Proof of the correlation identity for the honeycomb lattice 

From equation (|T2"j) 

< F{S)S, >= (F(S) HitfcoahipJijV) + SjsinhifUijV) + 1 - S?))f(x)\ x=0 (B.l) 



where, 



2_e^sinh{x) 

J[X) 2e^cosh{x) + V [a - Z) 



we obtain < S^Si >, for the honeycomb lattice, 



< S Q Si >=< Si(l + Si sinh JV + S?[cosh JV - 1]) 
x (1 + ^ sinh JV + Sf [cosh JV - 1] ) 
x (1 + 5 3 sinh JV + Sf [cosh JV - 1] ) > 

(B.3) 
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where Si, £2 and S3 are the neighbours of Sq. 
Or, 



< S0S1 >— 3ai < S\Si > +6(a 2 — ai) < SiS% > 

+a 3 < SiSaSs > +(ai - 2a 2 + a 4 ) < S'lSfSf >, (B.4) 



where, 



01 =sinhJV-/(x) \ x=0 =f(/3J) 
a 2 = sinh JVcosh JV • f{x) \ x=0 = 1/2/(2/3 J) 
a 3 = sinh J 3 V • f{x) | x=0 = l/4[/(3/3J) - 3/(/3J)] 
a 4 = sinh JV cosh 2 JV • /(z) | x=0 = l/4[/(3/3J) + f(/3J)] (B.5) 
From those results we obtain equations (fT!))) and (|17p of section 13.11 

Appendix C. Spin Correlation for the One-Dimensional S=l Blume-Capel 
Model 

For the linear chain, we have, 

< S >=< (1 + Si sinh JV + 5 2 [cosh JV - 1]) 
(1 + 5_i sinh JV + 5 2 ! [cosh JV-1]) > -f(x) | x=0 

(C.l) 

with /(x) given by expression (fTCTj) and Si and S-i are neighbors of So- We obtain for 
the two-spin correlation function 

< S S R >=< (SiS s + S-!S R ) > f(k) 

+ < (SiS-iSiSn + SiS^S^Sn) > [\f{2k) - f(k)) (C.2) 

where k = f3J. Applying the inequalities ( [l2| . [l3jj ) 

< SfS—iSfi S—iSr > 

< SI^Sr ><< 5i5 fl > (C.3) 

we get, 

< S Sr >< (< S 1 S R > + < S^Sr >)f(k) 

+(< S^Sr > + < SiSr >)[l/2/(2fc) - f(k)} (C.4) 

Defining C(R) =< S S R > we get 

C(R) = A(k)(C(R + l) + C(R-l)), (C.5) 

where A(k) = /(2fc)/2. 

If 7(i?) = C(R + l)/C(R) is inserted in the previous equation we get 

l = A(fc)( 7 (i?)+ 7 (i?)- 1 ). (C.6) 

So, C(R) = -f R and 

1 + y/l- 2/(2/3.7) 

7 = 7W ' (CJ) 

12 



Acknowledgements 

FCSB is grateful to CAPES/Brazil for the financial support that made possible his 
visit to the UFSJ/Brazil. ALM acknowledges financial support from CNPq-Brazil and 
FAPEMIG-Brazil. 



References 

Capel H W 1966 Physica 32 966 
Blume M 1966 Phys.Rev. 141 517 
Callen H B 1963 Phys.Lett. 4 161 
Fisher M 1967 Phys.Rev. 162 480 
Simon B 1980 Commun.Math Phys. 77 111 

Brydges D, Frohlich J and Spencer T 1982 Commun.Math.Phys. 83 123 
Lieb E 1980 Commun.Math.Phys. 77 127 
Siqueira A F and Fittipaldi I P 1986 Physica A 138 592 
Griffiths R B 1969 J.Math.Phys. 10 1559 
Sylvester G S 1976 J.Stat.Phys. 15 327 

Fernandez R., Frohlich J. and Sokal A. D. 1992 Random Walks, Critical Phenomena and Triviality 
in Quantum Field Theory (Springer- Verlag, Berlim) 
Braga G A, Ferreira S J and Sa Barreto F C 1993 Braz.J.Phys. 23 343 
Braga G A, Ferreira S J and Sa Barreto F C 1994 J.Stat.Phys. 76 819 
Szasz D 1978 J.Stat.Phys. 19 453 

Newman C 1975 Zcitschriftfur Wahrschcinlichkcits Theorie 33 75 
Ginibre J 1969 Phys.Rev.Lctt. 23 828 
Yukscl Y, Akinsi and Polat H 2009 Phys.Scr. 79 1 
Wilson K 1974 Phys. Rev. D 10 2445-59 

Seiler E 1982 Gauge Theories as a Problem of Constructive Quantum Field Theory and Statistical 
Mechanics, Lect. Notes in Physics 159 
Kogut J 1979 Rev. Mod. Phys. 51 659 

Sa Barreto F C and O'Carroll M L 1983 J.Phys.A 16 L431 
Glimm J and Jaffe A 1981 Quantum Physics (New York: Springer) 
Sa Barreto F C and O'Carroll M L 1983 J.Phys.A 16 1035 
Tomboulis E, Ukawa A and Windey P 1981 Nucl. Phys. B 180 [FS2] 294-300 
Griffiths R B 1967 J.Math.Phys. 8 478 
Griffiths R B 1967 J.Math.Phys. 8 484 
Griffiths R B 1967 Commun.Math.Phys. 9 121 

Plascak J A ; Caparica A A; Silva C J 2006 Phys. Rev. E 73 36702 
Xavier J C, Alcaraz F C, Pena Lara D and Plascak J A 1998 Phys. Rev. B 57 11575 
Beale P D 1986 Phys. Rev. B 33 1717 



13 



